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ABSTRACT 

Sensitive surveys of the Cosmic Microwave Background (CMB) will detect thousands of galaxy 
clusters via the Sunyaev-Zel'dovich (SZ) effect. Two SZ observables, the central or maximum and 
integrated Comptonization parameters y ma x and Y, relate in a simple way to the total cluster mass, 
which allow the construction of mass functions (MFs) that can be used to estimate cosmological 
parameters such as the ratio of the average matter density to the critical density Om> the normalization 
of the spectrum of initial density perturbations <7 8 , and the dark energy equation of state parameter w. 
However, clusters form from the mergers of smaller structures, events that can disrupt the equilibrium 
of intracluster gas upon which SZ-M relations rely. From a set of N-body/hydrodynamical simulations 
of binary cluster mergers, we calculate the evolution of Y and y ma x over the course of merger events 
and find that both parameters are transiently " boosted," primarily d uring the first core passage. We 
then use a semi- analytic technique developed bv lRandall et al.l (|2002l ) to estimate the effect of merger 
boosts on the distribution functions YF and yF of Y and y ma x, respectively, via cluster merger histories 
determined from extended Press-Schechter (PS) merger trees. The scatter in the Y—M and y m ax-M 
relations from merger boosts are found to be ~ 2% and 25-30% respectively. To determine Qm, &8, 
and w, the boosted and nonboosted YFs and yFs are fit with analytic PS distributions as a function of 
rcdshift. We find that boosts do not induce an overall systematic effect on YFs, and the values of Qm, 
as, and w (assumed constant) were returned to within 2% of values expected from the nonboosted 
YFs. The boosted yFs are significantly biased, however, causing VLm to be underestimated by 15-45%, 
erg to be overestimated by 10-25%, and w to be pushed to more negative values by 25-45%. We also 
fit YF as a function of redshift to cosmological models in which the dark energy parameter w varied 
with redshift to assess the effects of mergers on the inferred change in w with redshift. The values of 
CIm, &8, and the low-redshift value of w (u)q) were again reproduced to within 2%. For the largest 
change in w with z, which occurred between z = and z = 1 for the models assumed, it was increased 
by about 0.04. Although this is twice as large as the merger effect on a constant value of w, it is still 
reasonably modest. We confirm that the integrate d SZ effect, Y, is far more robust to mergers than 
Ymax, as previously reported by iMotl et al.1 (j2005) and similarly found for the X-ray equivalent Yx 
( Kravtsov et al.l l2006t IPoole et al.ll2007| ). and we conclude that Y is the superior choice for a mass 
proxy when using SZ observations of galaxy clusters to constrain cosmological parameters. 
Subject headings: cosmic microwave background — cosmological parameters — galaxies: clusters: 
general — hydrodynamics — intergalactic medium — large-scale structure of uni- 
verse 



1. INTRODUCTION 

The evolution of galaxy cluster abundance traces the 
massive end of the spectrum of initial density fluctua- 
tions and therefore is sensitive to cosmological parame- 
ters such as the ratio of the average matter density to the 
critical density Q.^ = 8ttG~p/ (3-ffg), the normalization of 
the power spectrum of initial density fluctuations as, and 
the dark energy equation of state parameter w, equal to 
the ratio of the pressure to the energy density of dark en- 
ergy. Here, Hq is the Hubble constant and p is the aver- 
age density in the universe. This sensitivity exists due to 
an exponential turnover at high masses in the mass func- 



tion (MF) of clusters, which can be pre dicted from a well- 
established theoretical framework (e.g. | Henrv fc Arnaudl 
ll99lHKitavama fc Sutdll996tlHaiman et al.ll2001f) . How- 
ever, only gravitational lensing, which remains observa- 
tionally challenging, directly measures the total mass of 
clusters. In order to get masses for the large number of 
clusters needed to construct the MF, it is often neces- 
sary to use a more observationally accessible quantity, 
such as the temperature or luminosity of X-ray emit- 
ting intracluster gas, from which the mass can be de- 
termined via some physical model. Relations between 
cluster mass and such a proxy typically require the gas 
to be in virial equilibrium; however, many processes are 
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known that can disrupt the gas, inclu ding cluster mergers 



(iRicker fe Sarazinl 



120021: iPoole etafl 



2001 



hereafter RS: lRitchie fc Thomas! 
2006)) and AGN jet-blown radio bub- 
20051) . 

many ongoing and planned surveys 
using the Sunyaev- Z erdovich (SZ) ef- 



bles (jMcNamara et al 

There are 
of cl usters 

feet dSunvaev fc Zel'dovichl fl97l iBirkinshawi Il999l) . 
which has the advantage of being effectively redshift- 
independent. The SZ effect is proportional to the integral 
of the electron pressure along the line of sight and can 
be characterized by the Comptonization parameter 
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n e T e dl cx / P„dl 
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where n e is the electron number density, T e is the elec- 
tron temperature, P e is the electron pressure, and I is 
the distance along the line of sight. The actual SZ flux, 
measured as a decrement or increment in the Cosmic Mi- 
crowave Background (CMB), depends on frequency and 
is subject to relativis tic effects for hi gh temperature plas- 
mas (for a review see lRephae li f 995). Because we do not 
want to tie our results to any particular observational 
project, we use the frequency-independent Comptoniza- 
tion parameter in the following study, as has been stan- 
dard in the literature. Also, we ignore any relativistic 
corrections as they are only relevant for the most mas- 
sive clusters and because they modify y in a complicated 
way that depends on frequency. 

In general, SZ observations will give an image of the SZ 
effect or y across the cluster. While specific values of y, 
for example the central or maximum value for a cluster 
(hereafter y ma x), are not expected to be a particularly 
robust proxy for the mass, the integrated Comptoniza- 
tion parameter Y disp lays a tighter correlation with mass 
(jReid fe Spergelll2006l ). This is defined as 
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where A is the projected surface area of the cluster on 
the sky, f2 is the solid angle, c?a is the angular diame- 
ter distance to the cluster, and the integral is over the 
entire cluster on the sky. Because the integrated Comp- 
tonization parameter is a global quantity, proportional 
to J P e dV or the thermal energy content of the elec- 
trons, it should be less sensitive to non-equilibrium pro- 
cesses, which tend to be more localized in cluster cores. 
The usefulness of SZ surveys to constrain cosmological 
parameters has a lread y been discusse d ext ensively (e.g., 
Carlstrom et al.ll200aTHaiman et all 120011: [Holder et all 
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As with X-ray proxies for mass, the regularity of an 
SZ-M correlation relies on the fact that many clusters 
are energetically close to equilibrium. However, dynami- 
cally unrelaxed clusters should add scatter to this corre- 
lation. One mechanism known to disrupt the gas is clus- 
ter mergers, a direct consequence of hierarchical struc- 
ture formation. How mergers affect the state of the gas 
will depend on the details of the individual mergers and 
their frequency, both of which depend on the cosmolog- 
ical model. To assess the utility of a mass proxy, such 
as the SZ effect, we need to quantify how mergers will 
affect the observed MF and consequently the estimation 
of cosmological parameters. 



Current cosmological simulations, which accurately 
trace the collapse of structure and thus the merger his- 
tory of clusters, cannot yet build the large samples of 
clusters at sufficient numerical resolution to constrain 
fundamental parameters and assess any potential bias 
due t o mergers - though t his approach is becoming viable 
(e.g.. lHallman et aT1l2007f h Typically, N-body cosmologi- 
cal simulations of dark matter are re-simulated to include 
various types of gas processes such as "preheating," ra- 
diative cooling, and AGN feedback, from which the scat- 
ter to an observed SZ-M correlation can be estimated. 
Depending on the resolution of the re-simulated hy drody- 
nami c grid, these studies p roduce samp les of ~ 10 ()Nagail 
1 20061 iBonaldi et all 120071) to ~ 100 (|Motl et al J 120051 : 
Ida Silva et al.|[2004|) clusters. Based on similar samples 
of simulations. iKravtsov et all (|2006l ) have defined an SZ- 
like X-ray observable, Yk, which they have shown to be 
robust to nonequilibrium gas physics with cosmological 
simulations. Though suited to understanding the physi- 
cal processes that add statistical scatter to SZ-M or sim- 
ilar relations, these samples are too small to assess the 
effect of the scatter on the determination of cosmolog- 
ical parameters, especially the effect of relatively rare, 
major merger events on the mass estimate of similarly 
rare massive clusters. To include these rare events and 
focus expressly on the role of mergers on SZ-M relations 
and cosmological parameter estimates, we take a semi- 
analytic approach that avoids simulating every possible 
merger within a cosmological framework. 

Specifically, we carefully examine the evolution of the 
SZ observables Y and y max for a discrete set of detailed 
N-body/hydrodynamical simulations of binary cluster 
mergers, generalize the results by identifying and pa- 
rameterizing the major transient features, or boosts, and 
then apply these boosts to the merger histories of many 
clusters generated semi-analytically via computationally 
c heaper merger trees. W e closely follow the methodology 
of Randal fet al.l (|200a hereafter RSR) , who similarly in- 
vestigated the effect of merger boosts on the X-ray ob- 
servable mass proxies Lx and Tx, the X-ray luminosity 
and temperature respectively, and the bias such boosts 
induce upon estimates of Qm and as from the inferred 
MFs. 

To assess the impact of a particular world model or 
cosmology on our results, we consider a "flat" cosmol- 
ogy with a cosmological constant, i.e. the ACDM con- 
cordance model, along with an "open" and Einstein-de 
Sitter ("EdS") world model for comparison; the relevant 
parameters are summarized in Table [1] The dark energy 
equation of state and its evolution are only examined for 
the flat universe. The Hubble constant is parameterized 
as 100 h km/s/Mpc throughout. 

In this paper, we assess the transient boosting of the SZ 
observables Y and 2/ max during cluster mergers, and the 
systematic influence of mergers on cosmological param- 
eter values derived from inferred cluster MFs. In §[2]we 
describe the binary cluster merger simulations from lRSl 
and the evolution of Y and y max during mergers. In § [3] 
we discuss the generation of cluster merger histories from 
merger tre es created via the extended Press-Schechte r 
formalism (|Press fe Schechterlll974]:lLacev fc Coldll993D . 
fit analytic functions that describe the transient behavior 
of merger boosts in the simulations, and generalize these 
functions to the entire family of possible mergers. In § [?] 
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the effect of boosts on the SZ-M relations are analyzed, 
in § [5] the distribution function proxies for the MF and 
the effect of boosts on them are described, and in § [5] 
the distribution functions are used to assess the impact 
of mergers on the cosmological parameters Qm, and 
w. Our results are discussed and summarized in § [7] 

2. MERGER SIMULATIONS 

To infer the effect of mergers on the SZ properties of 
clusters, detailed N-body/hydrodynamical simulations of 
every conceivable combination of cluster mass and im- 
pact parameter would be ideal. A realistic alternative is 
to u se a small but representative set of simulated mergers 
(|RSl ) and interpolate or extrapolate from them the ex- 
pected behavior of SZ observables for any set of merger 
parameters. 

A de tailed description of the simulations can be found 
in IRSI . Eight simulated binary cluster mergers were 
available with 3 mass ratios M > /M < = 1, 3, and 6.5 
each for 3 impact parameters b = (0, 2, 5)r s except the 
M > /M < — 6.5, b — 2r s case. Here r s is the scale ra- 
dius in the NFW p rofile for the more massive cluster 
(N avarro et al.l ll997). In all simulations, the less mas- 
sive cluster's mass was fixed at M < = 2 x 10 14 M©. Note 
that the M > jM < = 6 .5 simulation runs are not specifi- 
cally mentioned in IRSI . although they were generated by 
the same means as the other simulations. 

2.1. Equilibrium Y-M and y m:ix -M Relations 

To compare the SZ properties of merging clusters with 
those of similar clusters that are not undergoing merg- 
ers, we need an equilibrium SZ -M relation. The theo- 
retical models of clusters used in lRSl are designed to rep- 
resent observed, non-cooling flow clusters and to have X- 
ray temperatures typical of present day "rich" clusters. 
These initial conditions therefore include "preheating" 
and radiative cooling, though radiative cooling is ignored 
as a dynamic process as the cooling timescale is designed 
to e xceed a Hubble time. Though cooling is absent in the 
I RSI mergers, our results for y max generally agree w ith a 
similar set of cluster simulations (jPoole et al. 2007]) that 
do include radiative cooling. In any case, we are inter- 
ested in the change of Y or y max due to mergers and not 
in precisely characterizing the equilibrium state of clus- 
ters. To accurately assess the relative effect of m ergers 
on the SZ effect, we take the initial clusters in I RSI as our 
equilibrium clusters, which should correspond well to ac- 
tual clusters since they were built to resemble observed, 
relaxed clusters. 

Self-s imilar scaling relations derived from virial argu- 
ments (|da Silva et al.ll2004t [Cohn k Kadotal I2005D give 
Y tx M 5 / 3 f g , where M is the virial mass and f g is the 
gas mass fraction. For masses M > 10 14 M Q , f g oc M 1 / 3, 
though f g steepens at smaller masses. This general trend 
of increasing f g with mass has been observ ed for re- 
laxed, nearby clusters (Vikhlini n et al J 12006). Assum- 
ing all clusters have similar density profiles, we find that 
Umax Mfg. For the initial clusters in the simulations, 
we calculate exact solutions for Y and y max : 

„ m „ = 8 . 8 4xl0»( I5 4-) 2 (^)"/ s . (4) 



Here, r s and f g are found numerically (equations (20)- 
(23). lRSl ). Over the range of cluster masses we consider, 

Y and y m ax scale approximately as Y oc M 2 and y max oc 
M 1 - 3 . 

In practice, we fit the numerical solutions for Y(M) 
and ?/ max (M) each to a power law times a 13 degree poly- 
nomial. The high order of the polynomial is required pri- 
marily because we need the derivatives of the function to 
compute the Y and y max distribution functions (YF and 
yF respectively). The fractional error in the derivatives 
of the fits is < 1% for both Y(M) and y max (M), and 
better than that for the fits themselves. 

2.2. Merger Boosts to Y and ?/ max 
2.2.1. Generating Y and 7v max from the Simulations 

For each simulation in IRSI . the behavior of t he X -ray 
temperature and luminosity was calculated (see lRSl . Fig- 
ures 5 & 8) as a function of time. We would like similar 
curves for Y and y max ; however, these quantities were not 
calculated during the simulations, so we need to evaluate 
them from saved 3D "snapshots" of the simulation grid 
in order to recreate the evolution with time. For most of 
the runs, 40 to 60 snapshots were saved fairly regularly 
over the 14 Gyr the mergers were followed. From the 
gas pressure distribution, the Comptonization parame- 
ters can be calculated individually for each snapshot and 
combined to trace the evolution of Y and y max during 
the merger. 

Simulated SZ images for any orientation can be gen- 
erated for each snapshot. As an example, Figure [T] 
shows 100 x 100 pixel images from 2 snapshots of the 
M > /M < — 3, b — 2r s merger. For both of these images, 
our line-of-sight is oriented at 45° to the merger axis and 
rotated 45° azimuthally from the merger plane. In this 
particular example, the clusters are seen just before and 
just after the first core crossing, which generally corre- 
sponds to the maximum transient enhancement of both 

Y and y max . Note that while the images look qualita- 
tively similar, the scale of the image after core passage 
is twice that of the pre-core passage snapshot, suggest- 
ing that both Y and j/ max should get "boosted" during a 
merger event. 

To compute Y = J yd A = J P e dV, we simply 
add up the pressure in each computational cell weighted 
by the cell volume so that Y = ^ T c2 P^AVi, where 

i 

the sum is over all the cells in the 3D grid. We do not 
restrict the integration to the virial radius 7-200 or ^500 as 
in other cases where Y has been modeled ()da Silva et al.1 
I2004t iMotl et all 120051) for several reasons. First, the 
initial conditions for the simulated clusters cut off the 
pressure and density profile at the virial radius, so these 
definitions are at least initially equivalent. Also, during 
the merger there is no such well-defined radius as the 
gas is interacting violently. However, nearly all of the 
contribution to Y comes from gas inside the virial radius: 
99.5% initially and 95% after the clusters have merged 
and equilibrated. 

For each snapshot of each merger simulation, values 
of y m a X are computed for 339 orientations of the merger 
relative to our line of sight. Because the effects of the 
merger on the value of y max tend to vary the most near 
the merger axis, we more finely sampled the viewing an- 



4 



gles in this direction. The orientations sampled with re- 
spect to the merger axis are uniformly spaced in sin0, 
where 9 is the polar angle, such that Asin# = 1/15. The 
sampling of the azimuthal angle <f> is varied, to ensure 
relatively even spacing, as A0 = 8°/sin#. To determine 
J/max for each merger, snapshot, and orientation, values 
of y were computed by integrating along 16 lines of sight 
(equation ([1])), on a 4x4 grid, to form an SZ image of the 
cluster as seen from that orientation. The grid was then 
recentered on the maximum value of y and reduced in 
scale by a factor of 3.5, and y was calculated again. This 
procedure was repeated until the maximum value on the 
grid varied by less than 0.1% compared to the value from 
the previous iteration, and this y is adopted as y max . 

2.2.2. Correcting for Mass Loss Outside the Grid 

During each merger, some gas is flung out to large radii 
and lost from the simulation due to the finite size of the 
computational grid and outflow boundary conditions at 
the grid edge. Of course, once the gas is outside the sim- 
ulation grid, it is permanently lost. Noticeable amounts 
of gas do not leave the grid until after the first core pas- 
sage. Since we are mainly interested in the times when 
the merger boost is large, which occurs near the peak as- 
sociated with first core passage, our results are not par- 
ticularly affected by the lost gas. At late times, however, 
after the clusters have merged, Y remains below the ex- 
pected value for a cluster with mass M to tai = M< + M > . 
Since at these cluster masses Y cx M 2 and y max oc M 13 , 
we correct for the lost gas by taking 
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M gas (t = 0) 
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M gas (t = 0) 
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Here, Y" ca ic and y ca ic are the integrated and maximum 
SZ parameters calculated by integration over the grid 
prior to this correction. The correction is small; over 
the duration of the first peak in Y, which is much longer 
than the peak in y max , less than 5% of the gas has been 
lost from the grid. 

In fact, some of the lost gas exits the grid near to or 
above escape velocity, assuming a collisionless ballistic 
trajectory, so correcting for its loss may seem inappro- 
priate. The majority of the gas, except during the short 
period after the first core passage, effectively leaks out 
of the grid due to a lack of pressure support at the sim- 
ulation boundary. This artificially lower pressure prop- 
agates inward, requiring the correction we apply; other- 
wise, the boost will be slightly underestimated. After 
clusters have formed in cosmological numerical simula- 
tions, the gas fraction at the virial r adius is generally 
10% below the c osmic baryon fraction (|Crain et al.ll2007h 
lEke et aT] 1 19981 ). perhaps indicating that up to ~5% of 
the gas has been ejected, given tha t 5% of the baryons 
are in stars. The simulations of iRSl we utilize cannot ac- 
curately follow the merger to its true final state and so 
we cannot address the question of true gas ejection from 
clusters after merger events. However, the initial simu- 
lated clusters are constructed to match observed clusters 
with realistic gas fractions, so if gas is in fact lost, that 



effect is intrinsically included by IRSl and the resultant 
boosts in Y and y max . Our conclusions are not drawn 
from any late time evolution in the simulations, nor do 
we investigate the true post-merger state of clusters. 

Additionally, we correct Y for the slight evolution at 
large radii in the relaxed, pre-merger profiles of the sim- 
ulated clusters. Because the integrated Comptonization 
parameter is inversely proportional to a low power of 
cluster radius due to Y oc n e , the outer parts of a cluster 
contribute significantly to its overall value, as compared 
to Lx, which is proportional to rr e . The lower pressure 
in the outer regions can affect Y because there is more 
volume at large radii, even though ?/ max remains unaf- 
fected. We observe a slight drop in the pressure profile 
outside the central core over time before the individual 
clusters begin to interact, which is likely due to the arti- 
ficial truncation of gas at the virial radius - gas at this 
boundary is not in hydrostatic equilibrium in the sim- 
ulations and will flow outward, and the loss of pressure 
support will travel inward, readjusting the profile as the 
system tries to establish hydrostatic balance. While for 
the least massive cluster this effect is hardly noticeable, 
the magnitude of the effect increases with total cluster 
mass. Fortunately, the effect on Y(t) appears to be lin- 
ear in time, so we correct the time evolution of Y such 
that Y is forced to be constant before the clusters begin 
to interact, normalized to Y(t = 0). 

2.3. Evolution of y max and Y During Mergers 

In Figure[21 Y and y max are shown as a function of time 
for the merger simulations including the corrections de- 
scribed in § 12.2.21 For y max , the plot is shown for a 
viewing orientation at 90° to the merger axis and in the 
merger plane. The maximum boost for the head-on col- 
lision in y max is nearly a factor of 10, while the boost 
in Y is always less than a factor of 2, though the du- 
ration of_the_boost_^n Y is much longer than that for 
y max . iMotl et all (|2005| ) report a maximum boost factor 
in y max of 20 in cosmological simulations re-simulated to 
include gas hydrodynamics, twice the amount of boost- 
ing we find, though their result could be due to an ar- 
tificially high central temp erature in their pre-merger 
clusters (jLoken et al J [2002) . However, it is more likely 
the enhanced boost is due to the natural inclusion of 
multiple mergers and constant accretion along filaments, 
which are not included in binary merger simulations. 
For example, a triple merger between 2 equal mass clus- 
ters and a third subcluster with a tenth of one of their 
masses should yield a boost factor of 20, extrapolating 
our results to such a case. Additional pressure due to 
bulk motions within the pre-merger clusters, producing 
stronger shocks, may also lead to a larger boost. Glob- 
ally, the temperature profile of the initial clusters in IRSl 
agrees well with those clusters a ssembled in cosmological 
simulations (jLoken et al.ll2002[ ). so the precise origin of 
the discrepency is unclear. However, our boost factors 
are confirmed in a recent s et of binary cluste r mergers 
(|Poole et al.ll2006h . in which iPoole et all <|2007T l find y max 
to be boosted by a factor of ~ 10 (see their Figure 7). 

Essentially, y max traces the densest parts of clusters, 
which are the cores. These remain reasonably intact until 
near the time of first core crossing, which makes the peak 
in 2/m ax relatively narrow. On the other hand, Y involves 
a sum of all the gas, so it begins to get boosted as soon 
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as gas at large radii starts to interact, long before the 
cores approach, and the boost lasts longer, as gas in the 
outer regions needs more time to re-equilibrate. The time 
evo lution of y max is qualitatively similar to that found by 
RSl for the X-ray temperature and luminosity, quantities 
that are also dominated by the cores of clusters due to 
the fact that the X-ray emissivity depends on the square 
of the density. 

The plots of the evolution of Y already indicate that 
this parameter will not be strongly affected by mergers. 
First, the boosts in Y are smaller than in j/ ma x- Second, 
the boosts are not large compared to the equilibrium ef- 
fect of increasing the mass. Assuming Y oc M 2 , the boost 
factor B needed to exceed the final equilibrium value of 



Yfi 



is B > 



Y1+Y2 



M'(+M% 



For equal mass merg- 



ers, this condition gives B > 2, and from Figure [2] it is 
clear the boost factor is always < 2. If Y is used as a 
proxy to determine the mass of a cluster, the resulting 
value during the merger will nearly always lie between 
the individual initial masses of the subclusters and the 
final total mass. In a certain sense, this only affects the 
definition of when the cluster has merged, and the ap- 
plicable mass, and does not represent a real bias. We 
find that other mass ratios can boost Y beyond the fi- 
nal equilibrium value, but only by factors slightly larger 
than unity. 

3. MERGER TREES 

Structure formation and evolution are most easily 
traced through the mass function (MF) of dark matter 
halos, n(M, z), where n(M, z)dM gives the number of ha- 
los per unit comoving volume with masses in the range 
M — > M + dM. Currently, the MF for a given cosmology 
at a given redshift can be fou nd most accurately f rom nu- 
merical N-body simulations (Springcl ct al. 2005). While 
accessing the result s of these simulations ha s become 
more feasible (e.g., iLemson fc Springe! 12001 . a semi- 
analytic approach to obtaining the MF proves more prac- 
tical, especially since we are concerned with the relative 
effect of merger boosts on the underlying MF and not the 
precise nature of the MF itself. We follow the PS formal- 
ism, which agrees with the MF found in numerical sim- 
ulatio ns, especially at higher masses (|Brvan fc Normanl 
1998); specifically, we use extended P ress- Schechter the- 
ory as developed in lBond et al.l (|199fl) and lLacev fe Cole] 
(1993) and applied in |RSRj. Though the PS formalism 
fails to reproduce the MF found in numerical simula- 
tion s at very high redshifts and low cluster m asses (see, 
e.g. JSheth fc Tormedll999t iLukic et al.ll2007t ). it is more 
than sufficient over the redshifts (z = — > 2) and masses 
( M = 10 14 -> 10 16 M ? ) of interest here. 

iPress fc Schechteii(jl97l give the MF at some redshift 
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where a(M) is the current rms density fluctuation within 
a sphere of mean mass M, and S c (z) is the critical linear 
overdensity required for a region to collapse at redshift 
Z. The derivation of this expression assumes that halos 
grow from Gaussian density fluctuations that have larger 
amplitudes on smaller scales. Structure then forms hier- 
archically, with small halos collapsing first and merging 



to form larger halos. In this scenario, the highest mass 
halos, observed as clusters of galaxies, form most recently 
and should be most affected by merger processes at the 
present day. 

From this extended P S form alism, we follow the pro- 
cedure outlined in § 3 of iRSRl . in which a "merger tree" 
is generated for a present day cluster. The merger tree 
traces the merger and accretion history of a cluster of 
mass M back in time. For each time step, a progenitor 
cluster of mass M v i is chos en from a probability distri- 
bution (|Lacev fc Coldfl993l equation (2.25)), and since 
we only consider binary mergers, the mass of the other 
progenitor cluster is given by M p2 = M — M p \ . We will 
use the notation M > and M< for the larger and smal ler 
masses of the subclusters in each binary merger. RSR 111 
§ 3.1 discusses the disadvantages of dealing solely with 
binary mergers; however, our set of simulated mergers 
does not address more complex mergers, so we have no 
good way to derive a boost for them. Also, boosts are 
most dramatic for near equal mass mergers, and in such 
cases additional merger participants will likely be much 
less massive and have a negligible effect on the resultant 
boosts. However, one result of ignoring multiple mergers 
is that the merger tree-derived MF tends to overestimate 
the analytic PS MF for z > 0. The progenitor cluster 
with mass M P 2 is not taken from the PS distribution and 
is generally overestimated, so that the high mass end of 
the MF is overestimated at the expense of the very low 
mass end. Since we concern ourselves with the highest 
mass clusters, our resulting MFs will lie slightly above 
the analytic prediction, as illustrated in Figure which 
we must take into account when fitting MFs in § [6] 

A large number of merger trees was created with a 
broad span of initial cluster masses, and the distribu- 
tion of the initial masses was weighted so as to give the 
present day mass function. From ensembles of merger 
trees for the cosmologies of interest, we can find the MF 
at any redshift, and at any redshift we have each clus- 
ter's merger history, which can be used to determine the 
merger boost in some observable — in our case Y and 

?/max- 

Merger trees are a simple and computationally cheap 
way to simulate structure formation for a particular 
world model. But they are limited in that they only spec- 
ify progenitor cluster masses and discrete time intervals 
during which the mergers occur. All the dynamics and 
other details of a merger, however, are absent from EPS- 
derived merger trees. The information needed to connect 
the trees to our merger simulations is the masses of the 
clusters, the impact parameter b of the encounter, and 
the time of first core passage, which we designate as the 
time of the merger, t merge , in the merger trees. While the 
masses are provided by the merger trees, an appropriate 
b must be selected fo r each merger in the trees. We follow 
dM the method in § 6 of|RSRj, where a value for the spin pa- 
rameter is chosen from a Maxwell-Boltzmann-likc distri- 
bution, which represent s the observed dist ribution from 
numerical simulations (Bu llock et al.l |2~001[ ). all owing b 
to be derived from the chosen spin parameter (|Sarazinl 
2002). To determine the precise value of f mer ge, we sim- 
ply select a random time within the small discrete time 
step used in the merger trees, and take that time to be 
the instant of first core passage, since the merger could 
have occurred at any point within that time. 



6 



3.1. Merger Boost Histograms 

As discussed in iRSRl . the effect of a merger boost on 
a cluster whose history is characterized by a merger tree 
can be determined from a histogram which gives the mag- 
nitude of the boost as a function of time. Since the 
merger trees give a statistical description of the history 
of cluster mergers, it is sufficient to determine the distri- 
bution histogram of boosts versus the observed time i bs- 
The form of the histogram reduces the details contained 
in the curves in Figure [2] to a simpler, one-to-one func- 
tion that can be fit by the merger parameters M< , M> , 
and b. In the fits, we scale the impact parameter b by the 
core radii of the two merging clusters, b' = b/(r c< + r c> ) 
and the time by the ratio of the virial radius of the more 
massive cluster to the gas sound speed, t sc , A more de- 
tailed explanation of these scalings is given in § 5.3 of 

IrsM 

3.1.1. Fitting Y Histograms from Simulations 

The left panel of Figure [3] shows the cumulative time 
spent by the system above any given value of Y for the 
M > /M < = 1 merger simulations. We use cubic spline 
interpolation of the boost curves in the left panel of Fig- 
ure[5]to produce a smoothly varying histogram. In RSR, 
the Tx and Lx boost histograms were well-fit by hy- 
perbolas parameterized by the equations given in their 
Appendix B. We find hyperbolas also well-describe the 
Y his tograms, and we use the same parameterization as 
RSR with only a minor change given in Appendix [AJ 

3.1.2. Fitting j/ m ax Histograms from Simulations 

The procedure for ?/ max is slightly more complicated 
due to the orientation-dependence of the central Comp- 
tonization parameter. The evolution of the maximum 
value of y as a function of viewing angle varies more 
dramatically near the merger axis than perpendicular to 
it. As noted in § 12.2.11 for each merger simulation, the 
evolution of j/ m ax with time was calculated for 339 ori- 
entations, sampling more finely around the merger axis. 
Because y max really traces the cluster cores, the peaks 
in the curves are larger and have shorter durations. As 
a result, we found that simple interpolation did not suf- 
ficiently sample the peaks, so we use a superposition of 
Gaussians to fit the shape of the boost as a function of 
time. 

The merger trees contain no information about the 
orientation of the cluster mergers, and we assume an 
isotropic distribution relative to our line-of-sight. For 
our grid of 339 viewing angles, the probability of any 
one orientation is determined by the solid angle of that 
grid cell. We weight each orientation by this solid angle 
divided by An. All 339 y max {t) curves, weighted by their 
probability of being observed, are used to construct a his- 
togram like those described in § 13.1.11 The histograms 
for the M > /M < = 1 runs are shown in the right panel 
of Figure [3j Since these histograms include the distribu- 
tion of merger boosts for all orientations of the line of 
sight, the boosts need to be normalized to the pre-boost 
value of y max for some fixed orientation. The boosts in 
the right panel of Figure [3] were taken relative to the 
pre-boost y max as observed 90° to the merger axis and in 
the merger plane. Note that this is the same orientation 
assumed in the right panels of Figure [2] 



Because the y max histograms include the results from 
many different orientations, the high boost ends of the 
histograms decline more slowly with time than for Y or 
Lx or Tx- A different function is thus used to fit these 
histograms (see Appendix IA"|). 

3.2. Generalizing Merger Boosts for Arbitrary Mass 
Ratio and Impact Parameter 

As in lRSRl . the parameters of the fits to the boost his- 
tograms were fit to simple functions of the masses and 
impact parameter in the merger. The forms of these 
functions were chosen so as to have the correct asymp- 
totic forms (e.g., in the limit of large M>/M<). The 
free parameters of these functions were chosen to best fit 
the histograms from all 8 simulation runs. The values 
of these parameters are given in Table [5] (below in the 
Appendix) for the Y and y m ax histograms. 

The maximum fractional error in the fits to the boost 
simulation data for Y is < 3% except for the 2 runs 
with M > /M < = 6.5. Here, the evolution of the pressure 
distribution in the more massive cluster before collision 
dominates the time evolution of Y (see § 12.2. 2|) . The fits 
overestimate the boosts for the M > /M < = 6.5 simula- 
tions; however, the boosts themselves are small in this 
case, and the errors are still < 10%. 

For the j/ max fits, the average fractional error is typi- 
cally 4%, and the maximum error is < 10%. We found 
that the time sampling for the M> /M< =3,6 = 2r s sim- 
ulation run was too sparse around the boost to strongly 
constrain the shape of the j/ max histogram, so we did not 
use this run in our fits. 

3.3. Adding Boosts to Merger Trees 

With the fitted forms for the histograms for the 
strength of a boost versus time as a function of the masses 
of the merging subclusters M < and M > and the impact 
parameter b, the boosted values of Y or y max can easily be 
found for clusters from their past merger histories given 
by the merger trees. For any redshift or observed time, 
t bs, we search back through a cluster's merger tree and 
for every merger event, we find the boosted value of Y or 
Umax f° r that merger. If the boosted Y or y max exceeds 
the value given by our equilibrium equations and (H]) 
for the mass of the cluster at i bs, then we assign the 
boosted value to that cluster's observed Y or y max ; oth- 
erwise it acquires its equilibrium value. Boosted values 
less than those given by the equilibrium equations are not 
allowed because the analytic fits from which boost fac- 
tors are derived poorly describe the histograms, such as 
those shown in Figure [3J for negative and small positive 
boosts. While the discrepency between the simulation- 
based histograms and the analytic fits for small boosts 
leads to an underestimate of the number of these clus- 
ters, we are primarily concerned with the more dramatic 
effects caused by large boosts, which are well-described 
by the fits. 

4. SZ VERSUS MASS CORRELATION 

Once clusters observed at some redshift are assigned 
values of Y and y max based on each cluster's merger his- 
tory, we can evaluate the robustness of the Y — M and 
ymnx-M relations. The top panels of Figures [4] and [5] 
show Y and j/ max versus mass for clusters in our merger 
trees at z = and 2=1. Most clusters have nearly 
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unboosted values of Y and y max , while the number of 
clusters that deviate from either SZ-M relation drops 
roughly exponentially with the strength of the boost. 
We find that ~ 15% of clusters are boosted in y max by 
> 15% and in Y by > 0.1%. Note that the scatter in Fig- 
ures E] and [5] is due entirely to merger boosts and does 
not include observational error or scatter related to other 
physics. 

As expected, many clusters are found to have signif- 
icantly boosted values of y max , which overestimate the 
actual masses. However, there are almost no "boosts" 
to Y in Figured] Instead, we see clusters scattered be- 
low the Y—M relation, as is also seen, though to a lesser 
extent, in Figure [5] Clusters that fall below the Y-M 
relation were "observed" after a merger (after the peak 
of the boost), but before virialization. It should be noted 
that, according to ^ 13.31 a boost is only applied if it gives 
a Y or j/ max greater than its equilibrium value before the 
merger, while in Figures |4] and [5] the mass is taken to 
be the final, or merged, mass of the clusters. So, though 
clusters can never fall below their pre-merger equilib- 
rium relation in our formulation, a boosted cluster may 
fall below its post-merger value. The scatter below the 
Vmax-M relation is less pronounced due to the shorter pe- 
riod when the SZ effect is below the eventual equilibrium 
value (see Figure [2J . This feature is a general character- 
istic of observing a recent post-merger cluster and will be 
difficult to identify as such in an actual survey, and will 
likely affect the normalization of either SZ-M relation. 

In order to quantify the effect of mergers on the SZ 
versus mass relations, we fit power-law functions of the 
form 

or 

/ M \ a 

to all of the clusters with Y > 10~ 5 /i -2 Mpc 2 or y max > 
10~ 5 . We estimate the scatter with respect to the best 
fit, cifit, and also the scatter and offset with respect to 
the actual equilibrium relations for the SZ effect (equa- 
tions ([3]) and (TJI), °cq- We define the scatter as 

^2(Vi -y&t,i) 2 /y&t,i 



and <Teq is similarly defined, except N — 1 is replaced by 
N. The coefficients A and a of the fits along with the 
scatter are given in Tabled) The subscripts "b" and "nb" 
refer to clusters including merger boosts, and not includ- 
ing these boosts (where the SZ properties are given by 
the equilibrium relations). Note that we consider a log- 
arithmic distribution as in §[5] for the cluster masses, so 
that the fits in Table [H as well as the points in Figures 2] 
and [5l do not reflect the actual MF of clusters. 

Because the relative strength of boosts is mainly a 
function of mass ratio and is only weakly dependent 
on the absolute masses of the merging clusters, clusters 
are boosted somewhat uniformly in Y and y max across 
masses, which tends to change the normalization of the 
fit, A, and only to a lesser extent the slope, a. The inclu- 
sion of merger boosts, in the case of y max , could either 



flatten or steepen the slope. The local mass function 
is flatter at the low-mass end, so low-mass clusters ex- 
perience more high-mass-ratio mergers overall than high- 
mass clusters, thus flattening the y max -M relation. How- 
ever, when both minor and major mergers are consid- 
ered, at any given time the high-mass clusters are un- 
dergoing more merger events (see the relative change 
in yF over time for low- and high-mass clusters in Fig- 
ure [S]). Thus, at any given time a higher-mass cluster 
has a greater probability of finding itself in the midst of 
a merger of some type. If the mass function is oversam- 
pled at the high-mass end, this effect tends to steepen 
the y ma x-M relation. Because our cluster sample has a 
uniform distribution in log mass, we oversample the high- 
mass end relative to the low-mass end. Consequently 
our y max -M relation does not exhibit the flattening that 
we would expect if our cluster sample had been drawn 
from the correct mass function. As Figure \5\ shows, both 
high- and low-mass clusters exhibit the same number of 
large boosts, but the total number of boosted clusters is 
greater at higher masses. Most of the high-mass clusters 
with boosts have small boost factors that are difficult to 
see in the figure. 

The boosted normalization for Y is systematically 
lower than the nonboosted A, but by < 1%. In the case 
of y max , the normalization increases by ~ 10%. The off- 
sets to y max are due as much to clusters with small boost 
factors as to the rarer cases with very large boosts. Note 
that these clusters tend to be undergoing weaker mergers, 
which may be hard to detect. Thus, it may be difficult 
to expunge these clusters from SZ surveys, and the sys- 
tematic shift in the y max versus mass relation may bias 
cluster samples. The merger-induced scatter to the fit, 
CTflt, is ~ 2% for the Y-M relation and 25-30% for the 
Vm&x-M relation, and is nearly independent of the cos- 
mological world model and redshift. The scatter relative 
to the equilibrium relations, cr cq , increases with redshift 
since the merger rate is higher in the past, whereas the 
addition of boosted clusters adjusts the normalization A 
to minimize am, so the scatter remains about constant 
between redshifts. Also, because there are fewer clusters 
that show boosts in Y, the scatter <jgt is dominated by 
deviations of the equilibrium relation from a power-law 
form, which explains why a eq tends to be slightly smaller 
than ct^ for Y. 



5. DISTRIBUTION FUNCTIONS OF Y AND Ymax 

We computed the distribution functions for Y and 
J/m ax , which we refer to as the YF and yF, respectively. 
The distribution function YF is n(Y, z), where n(Y, z) dY 
gives the number of clusters per unit comoving volume 
at redshift z which have integrated SZ parameters in 
the range Y — > Y + dY. The yF distribution func- 
tion n(y max , z) is defined in an equivalent manner. To 
build the YF from a merger tree, we find all the clus- 
ters that exist at the "observed" redshift and assign a 
value of Y according to § 13.21 for the non-boosted YF 
and § 13.21 for the boosted YF. A cluster with integrated 
Comptonization parameter Yi is then added to a pre- 
determined bin YFj such that Y, bin < Y t < Y>j$ and 
appropriately weighted to convert the actual initial dis- 
tribution of z = cluster masses used in the merger trees, 



dN/dM , to the Press-Schechter distribution np$: 

n PS (M°,Z = 0) 



YF,- = YF, 



dN 



(yUn _ Y U a) 



(11) 



The initial distribution dN/dM° is logarithmically 
spaced to ensure good statistics at the high mass end, 
where clusters are rare, and to avoid creating an exces- 
sively large number of merger trees. 

For the non-boosted case, the YF or yF can be found 
directly from the equilibrium relations (equations ((3]) & 
10} ) and 



n PS (Y, z)dY = n PS (M, z)^dY 



n 



n PS (M,z) 



dM 



-dyn 



(12) 



(13) 



with nps(M, z) from equation ([7]). The derivatives are 
found from fits to the equilibrium relations (see § 12. ip . 

The agreement between the nonboosted merger tree- 
derived YFs (yFs) and the analytic Press-Schechter YFs 
(yFs) is shown in Figure [6] for the flat world model. Note 
that the merger trees seem to slightly overestimate the 
number of lower Y or y max (i.e. lower mass) clusters at 
higher redshifts, which is due to a feature of our merger 
tree p rocedure discussed in § [3] of this work and § 3.1 of 

Irsm 

The nonboosted and boosted YFs and yFs are also 
compared in Figure [6] The boosted YFs are almost iden- 
tical to the nonboosted YFs. The deviations from the 
nonboosted YFs are not systematic and are typically of 
a few percent and only visible in the residual plot. The 
boosted yFs, however, lie systematically above the non- 
boosted yFs at all 3 redshifts considered. The fractional 
deviation increases with both cluster mass and redshift. 
The increase with cluster mass shows that rare events 
involving major mergers of moderate mass clusters com- 
pete in frequency with the number of rare, very massive 
clusters with large equilibrium values of y max . The in- 
crease in the bias with redshift is apparently due to the 
higher merger rate in the past. Clearly, clusters with all 
values of z/ max are getting boosted to higher y max bins 
in the yFs over our considered range of y max , which in- 
cludes only the most massive clusters. Such a significant 
and systematic bias in the yF will affect estimates of cos- 
mological parameters, as discussed below in § [51 

6. DETERMINING COSMOLOGICAL PARAMETERS FROM 
THE MERGER TREE YFS AND YFS 

Although mergers strongly affect the SZ signals of a 
small fraction of clusters, because of the exponential 
high-mass drop-off in the YF and yF the effect of mergers 
on cosmological model fits to these distributions may be 
significant. To quantify this effect, we derive fits based on 
the analytic predictions of equations (fT2"j) and (flT?)) to the 
YF and yF using both boosted and nonboosted merger 
trees. The differences between the best-fit cosmological 
parameters derived in the two cases provide an estimate 
of the systematic bias introduced when merging effects 
are neglected. 

6.1. Varying only 51 m o-nd as 



For our 3 cosmological world models, we treat the 
binned YFs and yFs from the merger trees as observa- 
tional data and find best-fit values for the parameters VLm 
and as in equation (JT2J) or ([TBI . Due t o a ne ar degener- 
acy between fijvf and as (Bahcal l fc Fanlfl99 S?) at a single 
redshift, we simultaneously fit YFs and yFs for two red- 
shifts: at z = and at either z = 0.5 or z — 1.0. While 
in practice SZ surveys will observe clusters in a continu- 
ous range of redshifts, choosing only 2 redshifts simplifies 
the fitting procedure and illustrates the effect of merger 
boosts on these parameters. We choose to only fit clus- 
ters above a minimum value of Y min = 10 _5 /i _2 Mpc 2 
or y™ ax = 10~ 5 . These limits are consistent with the 
expected detection thresholds for upcoming SZ surveys, 
such as the AMI, ACT, and SPT projects fe.g.. lBartlettl 
|2006() . and the li kely confusion limi t for clusters with 
M < 1O 14 /i" 1 M (jHolder et al.l[20f)l . These limits also 
keep our fits from being biased by the large number of 
clusters at low masses. 

To evaluate the extent to which merger boosts affect 
the estimation of f2j\/ and as, we compare their fitted 
values from the boosted YFs and yFs to the fitted values 
from the nonboosted YFs and yFs. We do not com- 
pare best-fit parameters to the values used to create the 
merger trees because the trees tend to slightly overesti- 
mate the MF, an effect which increases with redshift and 
is discussed in § [Sj However, since we are only inter- 
ested in relative changes to the YF or yF due to boosts 
in Y or y max , this bias in the MF does not affect our 
results, though the best-fit parameters found from the 
nonboosted YFs or yFs may differ from the parameter 
values used to create the trees. Also, any bias caused 
by our chosen fitting method is accounted for by directly 
comparing the two YFs or yFs. 

The best-fit values of Q,m and a% for the flat, open, and 
EdS cosmological world models are given in Table [5] for 
both Y and y max . The parameter values used to create 
the merger trees are summarized in Table[T]for reference. 
In general, the results are independent of world model; 
cosmological parameter fits tend to be biased in the same 
direction by about the same amount. However, boosts 
to Y have almost no effect on fits to VLm and a$; the 
changes due to mergers are generally less than 1% and 
are not clearly systematic. 

In contrast, boosts to y max significantly bias the values 
of these parameters: Qm is underestimated by 15-30% 
and as is overestimated by 10-20%. The main effect of 
merger boosts is to increase the number of clusters de- 
tected in a particular yF^ bin; in other words, there is 
a systematic increase in the yF, as shown in the right 
panel of Figure [6] An overall increase in the normaliza- 
tion of the yF leads to an increase in the normalization 
of the spectrum of initial density perturbations, as- The 
total matter content, f^M, is also sensitive to the normal- 
ization, but it is near ly degenerate with as ~ 0.6f2 



-1/2 
M 



(|Bahcall fc Fan|[T998f) at a given single redshift. How- 
ever, fijyf is more sensitive to the change in the yF over 
time - the greater the density of matter, the faster struc- 
ture will grow. If various cosmologies with nearly identi- 
cal yFs at z — are considered, those cosmologies with 
smaller values of Q.m (and thus larger values of as) would 
produce yFs at z > that lie above the yFs of cosmolo- 
gies with larger 17 m values. As described in § [5j merger 
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boosts raise the yF most strongly at higher redshifts, 
so the change in the yF from one redshift to another 
is smaller than for nonboosted yFs, indicating a slower 
structure growth rate and therefore a smaller f2jw- The 
overall effect of mergers seems to vary with redshift; £Im 
and (is are found to be less biased when utilizing the yF 
at higher redshift (z = 1) even though this yF is frac- 
tionally more biased than the yFs at z = or z = 0.5. 

6.2. Fitting the Dark Energy Equation of State 
Parameter w 

Clusters of galaxies have been used to constrain the 
equation of state parameter w of dark energy, and there 
are extensive plans to improve th ese measurements in 
the future using SZ surveys (e.g., lHaiman et al.l 120011 : 
IWeller et al1l2002fl . In the ACDM flat world model, dark 
energy is assumed to take the form of a cosmological con- 
stant, which has a fixed w = — 1. Here, we assess the 
effect of mergers on the determination of w by allow- 
ing w to vary along with f2j\/ and as in fits to the flat 
world model YFs and yFs, following the same procedure 
outlined in § 16.11 We need new analytic, nonboosted 
Press-Schechter YFs and yFs that incorporate w ^ — 1, 
which we write as nps(Y, z, w) and nps(y max , z, w). The 
same basic form of nps can be generalized to a con- 
stant io/-l and a slowly varying parameterization of 
w(z) — wq + wia(l — a) = wq + w\zj (1 + z) 2 , where a 
is the scale factor and wq and w\ are constants and w\ 
is small. In a flat (£Im + &>de = 1) universe, we change 
the expression for the growth function D(z) as given in 
Appendix A of iRSRl and correct the comoving volume 
element dV such that 



/ dN \ w _ / dN y=-iMD(*,«)) /d|=-iy 

\dVdYJ ~ \dVdYj Wa / 



fdV_ 



dV u 



(14) 

where lMD[z ' w)) is n PS (Y, z) for the ACDM 

cosmology, but with the critical overdensity S c given by 
the new growt h function D(z,w) [equation (20) from 



the new growt h tunction L/(z,w) [equation (ZU) trom 
Percivall (120051) for consta nt w or equation (14) from 



Wang fc Steinhardtl (|1998l ) for for w(z)}, and d A is the 



angular diameter distance. The ratio of volumes is 



dV h 



dV" 



E(z, w) 



E(z,w= -1) 



(15) 



where E(z,w) = [fi M (l + zf + n DE (l + z )3+3»]i/2_ 
The same expression applies for the yFs by replacing Y 

(d w= ~ 1 \ 2 
-^p — ) from equa- 
tion pi)) . 

6.2.1. Constant w 7^ — 1 

When we allow for constant values of w that are not 
necessarily equal to —1, we find results qualitatively sim- 
ilar to what was found previously when only Qm and as 
were varied. Again, the boosted YFs give back nearly 
identical values for all 3 parameters to within < 1%. For 
2/max, merger boosts are found to bias the fitted values 
for Qm and as even more strongly, underestimating Qm 
by 30-45% and overestimating as by 20-25%. Also, w is 
found to be more negative in the boosted yFs by 25-45%, 
making y max a poor proxy if one aims to constrain the 



nature of dark energy. These results are summarized in 
Table H 

In the case of 2/ max , the boosted yFs favor more nega- 
tive values of w due to w;'s impact on structure formation. 
The yF is overestimated to a greater extent at larger red- 
shifts (see Figure |6|), which mimics more structure in the 
recent (z < 1) past. In turn, the appearance of more 
collapsed structures in the past relative to the present 
time implies that recent structure formation was slower 
than it actually has been, and that structure formation 
in the far past was correspondingly faster. In general, if 
we compare the effect of different values of w on structure 
formation by holding the present yF fixed, a more neg- 
ative w is better able to slow down cluster formation at 
later times as the strength of dark energy grows with the 
scale factor a since Vt DE = n DE ^ + z f {1+w) / E 2 (z). If 
cluster formation is slowing at the current epoch, when 
dark energy has recently become dominant, there must 
be more clusters in the recent past compared to the yF 
of clusters under the influence of a less negative w. 

A more negative w allows for even smaller values of 
f2j\/ to be fit to the boosted yFs, compared to its best-fit 
values when only Qm and as are varied. By anchoring 
the current yF, a more negative w decreases the influence 
of dark energy in the past, so Qm does not need to be as 
large to form the same amount of structure. The dark 
energy equation of state does not as directly affect the 
overall normalization of the yFs, so the bias to as remains 
consistent with the fixed w = — 1 fits. 

6.2.2. Slowly Varying w(z) 

If dark energy is not due to a cosmological constant, 
then it is possible that its equation of state might vary. 
"We have also determined the effect that merger boosts 
"cajn ' have on the SZ determination of the evolution of 
dark energy. We only consider the effect of boosted YFs 
in this section, due to the difficulty of using yFs to pin 
down even constant values of w. Choosing the parame- 
terization of w = wq + wiz/(l + z) 2 , where wq and Wi 
are constants, we determined 51m, wq, and w\ by fit- 
ting the boosted and nonboosted YFs. The validity of 
the form of the growth function we use for a flat universe 
requires that | 1 -C 1 _\i M , which implies that Wi C 1 
for w w -1 and H M w 0.3 (|Wang & Steinhardtl [l99l) . 
We do not constrain our best-fit value of W\ according to 
this requirement, however, nor do we consider any other 
parameterization of w. 

We found that w\ was not well-constrained by fitting 
the YFs or yFs at only two redshifts. Thus, we simul- 
taneously fit the distribution functions at the three red- 
shifts z — 0, 0.5, and 1. As in the constant w case, the 
boost-derived values of Am, ct§, and wq deviated from 
the nonboosted values only slightly, by +0.5%, —0.2%, 
and +2%, respectively. The best-fit values of w\ in- 
creased by 0.15 from the nonboosted value of —0.19 to 
a value for the boosted YF of —0.04. For the assumed 
variation of w with z, the largest change in w occurs be- 
tween the present time (z = 0) and z = 1; that change 
is Aw — wi/4. Thus, the merger boost effects on YF 
alter the maximum change in the w by about 0.04. This 
is about twice as large as the effect on wo, but is still 
relatively small. 

7. DISCUSSION AND SUMMARY 
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We have determined the effects of cluster mergers on 
their SZ properties, particularly the integrated Y and 
maximum j/ max Comptonization parameters. From a set 
of hydrodynamical/N-body simulations of cluster merg- 
ers, we determined the evolution of Y and y max over the 
period of interaction for mergers of various mass ratios 
and impact parameters, and we found that mergers tem- 
porarily "boost" both Y and y m ax- For y max , the boosts 
can be as large as an order of magnitude, although they 
occur for a short time (typically about half the sound 
crossing time of the cluster) , with the largest boosts oc- 
curring near the time of first core crossing. For major 
mergers, the boosts in the maximum Comptonization pa- 
rameter generally exceed the increase in j/ max when the 
systems have come into equilibrium. 

On the other hand, the boosts in Y are smaller (less 
than a factor of two), although they last longer (about 
two sound crossing times). Most importantly, the boosts 
in Y for major mergers are smaller than the increases 
in Y when the merged clusters have come into equilib- 
rium. Thus, one can think of the merger "boost" in Y 
as representing a stage in the evolution from two sep- 
arate equilibrium values to the final merged value, and 
not really being a "boost" at all. A simple physical ar- 
gument explains why the transient boosted values of Y 
are smaller than the final equilibrium values. From equa- 
tion ((21), it follows that Y is just proportional to the total 
thermal energy content of the electrons in the clusters, 
or just the total thermal energy if the electrons and ions 
are in equipartition. Now, a cluster merger involves the 
conversion of the bulk kinetic energy of the merging clus- 
ters into thermal energy. When the merger is complete, 
there is very little bulk kinetic energy remaining (per- 
haps, weak rotation or turbulence). Thus, one expects 
the thermal energy content of the merging clusters to 
be largest when they have achieved (or nearly achieved) 
equilibrium. Thus, the final equilibrium value of Y will 
tend to be larger than any transient value during the 
merger. 

We generalized the SZ boosts to mergers of arbitrary 
mass ratio and impact parameter and traced the merger, 
and thus boost, history of clusters with redshift using the 
EPS merger tree formalism. In general, merger boosts in- 
duced a relatively small scatter, ~ 2%, below the equilib- 
rium Y—M relation, while mergers induced a large scat- 
ter of 25-30% above the ym ax -M equilibrium relation. 
Power-law fits to Y and 2/ max as a function of mass show 
that while boosts do not affect the slope of the fit, the 
normalization was lowered by < 1% for Y and raised by 
~ 10% for 

We used the merger trees to derive the distribution 
functions of SZ parameters, YF and yF. We found that 
the boosted YF was not significantly biased relative to 
the nonboosted YF, while the boosted yF was strongly 
biased above the nonboosted yF for all redshifts. In gen- 
eral, the size of the merger-induced bias increased with 
redshift and with cluster mass. 

Using the YFs and yFs, we determined the best-fit 
values for the cosmological parameters Qm and cr% for 
the flat, open, and EdS world models, and also the dark 
energy equation of state parameter w for the flat uni- 
verse. Comparing the best-fit values of CIm and as for 
the nonboosted and boosted YF, no significant difference 
(< 1%) was observed. In contrast, the boosts to the yF 



decreased the best-fit value of Qm by 15-30% for the flat 
and open world models and increased the best-fit value 
of as by 10-20% for all world models. These results stem 
mainly from an overall increase in the yFs, which pushes 
as to larger values, and a greater increase in the boosted 
yF at higher redshifts relative to lower redshifts, which 
pushes £Im to smaller values. Allowing for a constant 
w ^ —1 in the flat world model, no systematic difference 
in fitted cosmological parameters was found between the 
two sets of YFs, though the merger- induced bias to Umi 
as, and w was exacerbated when using the yFs. We also 
considered a time-varying w(z) for the YFs, for which 
fijif , fTg , and wq were found to be consistent with the 
previous results for a constant w ^ — 1, though boosts 
increased the best-fit value of the dark energy evolution 
parameter w\ by about 0.15. The largest change in w 
occurs between z = 1 and z = in this model; thus, the 
change in w might be affected by as much as 0.04. This 
is about twice as large as the maximum change in the 
present-day value of Wo, but still is relatively moderate. 

These results agree with previous work which indi- 
cates that global observables such as Y or the equivalent 
X-ray/mass proxy Yx are more robust as mass proxies 
than the central or maximum Comptonization parame- 
ter. For example, from semi-analytic models of the intra- 
cluster medium (ICM) , iReid fc Spergell ((20061 ) generally 
find that Y tx f g M 5 ^ 3 , equivalent to our equilibrium def- 
inition of Y, with only a small scatter due to internal 
physics. A number of studies have used cosmological N- 
body simulations and re-simulated forming clusters with 
various kinds of g as physics to evaluate the scatte r in 
the y-M relations ((Nagail 120061 : iBonaldi et alJl2007f ). It 
is generally found that the normalization A varies signif- 
icantly depen ding o n the I CM physics, though the slope 
a does not. iNagail ((20061 ) reports a scatter of 10-15% 
in the Y-M re lation, much la r ger th an our scatter of 
- 2-3%. Also, iKravtsov etHl (|2006h defines an X-ray 
observable Yx — TxM &3S ^qq, which is similar to our 
Y oc J n e T e dV ~ T e J n e dV cx T e M gas ; they find a scat- 
ter in the relation of 5%-7%. 

While these studies intrinsically include mergers, they 
have limited statistics as they generally consist only of 
a small number of systems, ~ 10 or so. Some stud- 
ies have considered somewhat larger c luster simulation 
samples including h undreds of clusters ((Motl et al.ll2005t 
Ida Silva et al.l(2004l ) from cosmological simulations. Our 
results agree with their conclusions that the Y-M rela- 
tion is r elatively stable to mergers, unlike the y m ax-M 
relation. iMotl et all (|2005h find a scatter in their Y-M 
relation of 3-4% and in their y max -M relation of ~ 17% 
at z — due to mergers and other ICM physical pro- 
cesses. These results compare well with our scatter of 
2% and 24%, respectively. 

In a study similar to this work. lPoole et al. (2007) take 
a suite of binary cluster merger simulations to assess the 
effect of various observables, including SZ parameters, 
on scaling relations during mergers. The evolution of 
J/max (which they call y ) in their simulations is qual- 
itatively similar to our results in Figure [2] for various 
impact parameters and mass ratios. They also consider 
an integrated Comptonization parameter, but it is only 
integrated out to a radius 7-2500 and is thus much more 
dominated by core effects and not equivalent to our Y, 
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TABLE 1 

cosmological 
Parameter Values Used 
to Create Merger Trees 



Model flu Q\ as 



Flat 0.3 0.7 0.834 
Open 0.3 0.0 0.827 
EdS 1.0 0.0 0.514 



which is effectively integrated to at least r2oo, the virial 
radius. 

The large number of galaxy clusters expected from up- 
coming SZ surveys, both locally and at potentially high 
redshifts, heightens the prospects that clusters could play 
a decisive role in the era of precision cosmology, especially 
if the robustness of Y as a proxy for mass is confirmed 
in real cluster samples. 
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APPENDIX 
FITTING SIMULATION DATA 



We use the same basic forms and procedures to fit the merger boosts discussed in iRSRi . Appendix B. For the 
integrated Comptonization parameter, the bo osted part of the cumulative time distribution histograms is well-fit by 
hyperbolas similar in form to equation (Bl) of IRSRi with a slight modification: 
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TABLE 2 

Merger-Induced SZ-M Relations and Scatter 



Model 


z 


CToq 


Cfit 
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Y Flat 
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EdS 
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J/max Flat 
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0.292 
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9.01 
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24 
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24 
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9.46 






1.26 


EdS 
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0.290 


8.25 
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1 
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TABLE 3 

Best-Fit Values for Q m and <t 8 for Three World Models 



Model 




Boosts? | 


n M Di 


Terence 


cr 8 Di 


Terence 


Y Flat 


0,0.5 


no 


0.287 




0.857 








yes 


0.289 


0.7% 


0.854 


-0.4% 




0,1.0 


no 


0.277 




0.865 








yes 


0.277 


0.0% 


0.865 


0.0% 


Open 


0,0.5 


no 


0.278 




0.857 
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0.279 


0.4% 


0.856 


-0.1% 




0,1.0 


no 


0.279 




0.855 
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0.280 


0.4% 


0.855 


0.0% 


EdS 


0,0.5 


no 


0.932 




0.531 
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0.931 


-0.1% 


0.531 


0.0% 




0,1.0 


no 


0.874 




0.541 
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0.873 


-0.1% 


0.541 


0.0% 


Vmax Flat 


0,0.5 


no 


0.295 




0.844 
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0.199 


-33% 


1.020 


21% 




0,1.0 


no 
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12% 


Open 


0,0.5 


no 
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0.984 
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no 
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0.236 


-16% 
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13% 


EdS 


0,0.5 


no 


0.953 
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0.921 


-3.4% 


0.589 


12% 




0,1.0 


no 


0.905 




0.532 
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0.924 


2.1% 


0.590 


11% 



TABLE 4 

Best-Fit Flat World Models with Constant w 





Z 


Boosts? | 




Difference 


0"8 


Difference 


w 


Difference 


Y 


0,0.5 


no 


0.314 




0.837 




-0.879 






yes 


0.316 


0.6% 


0.835 


-0.2% 


-0.885 


0.7% 




0,1.0 


no 


0.275 




0.874 




-1.062 
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0.271 


-1.5% 


0.877 


0.3% 


-1.080 


1.7% 


Vmax 


0,0.5 


no 


0.324 




0.823 




-0.861 
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0.173 


-47% 


1.082 


24% 


-1.255 


46% 




0,1.0 


no 


0.279 
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-0.987 
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0.192 


-31% 


1.045 


21% 


-1.240 


26% 
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Fig. 1. — Images of the SZ parameter y from 3D snapshots of the 1:3 mass ratio, 2 r a impact parameter merger simulation. Here r s is 
the NFW scale radius of the more massive cluster. Both images are viewed from a line-of-sight which is rotated 45° from the merger axis 
and 45° azimuthally from the merger plane. Left: 386 Myr before first core crossing. Right: 114 Myr after first core crossing. 

Three parameters describe the function: the maximum boost [K/y(0)] poa k, the boost duration (t/t sc )Y, and the 
eccentricity of the hyperbola e. The fit values for these parameters between simulation runs co uld b e reproduced with 
the same functions of fractional mass increase /m and normalized impact parameter b' used in lRSRl . provided here for 
completeness: 
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+ 1, 



e(f M ,b') = 



C + b 12 



My)- H ln(M< /3 



+ Af* /3 ) 



I + b' 2 



(A2) 



(A3) 



(A4) 



As in the text, the impact parameter is scaled by the core radii of the two merging clusters, b' = b/(r c< + r c> ), M < 
and Af> are the masses of the less massive and more massive cluster (in M Q ), respecti vely, a nd the fractional mass 
increase /m = M < /(M < + M>). Motivations for these forms are given in Appendix B of RSR|. 

The variation of y max with the viewing angle of the merger causes the histograms of values of time versus j/ max 
to be broader than the histograms for Y (Figure [3]) . This difference makes hyperbolae a poor representation of the 
histogram shapes. We find a suitable replacement in another 3 parameter function 



h, i—)=P In (l - ^) - 

^sc / \ ypcak / 2 ?/pcak 



hi I ±. 

t-sc 



(A5) 



with similarly defined parameters for the maximum y max boost ^ ypca ^ , the power law slope P, and the boost duration 

"■(£),/ 

" ' A! " ' 1, (A6) 
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F + b' 1 - 5 
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ln(M< + M>) - H ln(A/< /3 + M> /J ) 
I + b' 2 



G 



(A7) 



(A8) 



The best-fit values found for A-I are given in Tabled] Note that A-I are found assuming that the value of y m ax(0) 
is taken along the merger axis, which is twice the value of y ma , x (0) used in Figures and [31 for which the value 
perpendicular to the merger axis is used. 
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Fig. 2. — Evolution of the SZ effect during a merger. In each panel, different curves are for different values of the merger impact parameter: 
6 = (solid, black line), 2 r s (dashed, red line), and 5 r a (dotted, blue line), where r s is the NFW scale radius of the more massive cluster. 
In the left panels, the b = 2r s simulation run is offset downward by 0.15 and the 5r s run is offset downward by 0.3 for clarity. The time 
is scaled by the sound crossing time t BC of the more massive premerger cluster. Left: Integrated Comptonization parameter Y versus time 
for the 1:1 (top), 1:3 (middle), and 1:6.5 (bottom) mass ratios. Right: Maximum Comptonization parameter y ma x versus time for the 1:1 
(top), 1:3 (middle), and 1:6.5 (bottom) mass ratios. The mergers are observed 90° to the merger axis and in the merger plane. 



TABLE 5 

Fitting Parameters for Merger Boost Histograms 



Boost ABCDEFGH I 



Y/Y(0) 95.69 0.8793 66.72 94.83 0.3621 173.3 33.36 0.2793 473.3 

3/max/ymax(0) 26.55 0.5776 4.052 6.310 - 4.569 2.250 1.785 13.76 
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YA(0) y™,/y™>(0) 

Fig. 3. — Left: Histogram of the total time the integrated Comptonization parameter Y is above some fraction of its initial premerger 
value Y(0), scaled by the sound crossing time t sc of the more massive premerger cluster. Histograms are shown for equal-mass mergers at 
3 impact parameters 6 = (solid, black line), 2 r s (dashed, red line), 5 r s (dotted, blue line), where r B is the NFW scale radius of the more 
massive cluster. Right: Histogram of times for j/max. 




1 10 1 10 

Mass (h~ 1 1O 15 M ) 

Fig. 4. — Integrated Comptonization parameter Y (top panels) versus total mass in the flat cosmology at z = (left panels) and z = 1 
(right panels) for clusters with Y > 10~ 5 h~ 2 Mpc 2 . The combined mass of both merging clusters is used if Y is boosted or t obs > t mer ge, 
where i mor ge is the time of maximum boost. The apparent solid line is the result of many individual clusters at or near their equilibrium 
values of Y. In the bottom panels, the ratio of the boosted clusters to their equilibrium values for each redshift is shown. Each panel 
contains 5190 clusters. 
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icr 5 icr 4 icr 3 icr 5 kt 4 icr 3 

Y (t-T 2 Mpc 2 ) y mo> 



Fig. 6. — Boosted (dashed line) versus nonboosted (solid line) integrated Comptonization parameter function YF (left panel) and 
maximum Comptonization parameter function yF (right panel) histograms for z = (top, black), z = 0.5 (middle, red), and z = 1 (bottom, 
blue) in the flat universe. The smooth curves are the analytic PS predictions at each redshift given by equation 101 • The residual plots 
give the difference in the logs between the boosted and nonboosted YFs and yFs. Note the significant difference in scales of the residuals 
between the YFs and yFs. 



